Brueckner-Hartree-Fock study of circular quantum dots 
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We calculate ground state energies in the Brueckner-Hartree-Fock theory for N electrons (with 
N < 20) confined to a circular quantum dot and in presence of a static magnetic field. Comparison 
with the predictions of Hartree-Fock, local-spin-density and exact configuration-interaction theories 
is made. We find that the correlations taken into account in Brueckner-Hartree-Fock calculations 
give an important contribution to the ground state energies, specially in strongly confined dots. In 
this high-density range, corresponding in practice to self- assembled quantum dots, the results of 
Brueckner-Hartree-Fock calculations are close to the exact values and better than those obtained in 
the local-spin-density approximation. 

PACS numbers: 73.21. La,71.10.-w 
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I. INTRODUCTION 

The Hartree and Hartree-Fock (HF) mean field ap- 
proaches have been extensively used in the past to 
study atomic-like properties of semiconductor quan- 
tum dots such as, e.g., those measured in conduc- 
tance and capacitance experiments^Si&l Symmetry re- 
stricted approximations, as well as spin and/or space 
unrestricted HF solutions have been analyzedASiZ&SiiS 
These methods give results that are satisfactory for 
a qualitative understanding of many properties of 
these systems. However, comparison with exact 
configuration- interaction diagonalizationiiii2iiiiiii5iiSiil 
(CI) and Quantum Monte Carlo 18 i 19 i 20 i 21 i 22 (QMC) stud- 
ies show discrepancies in the total energies that are sub- 
stantial on the relevant energy scale. 

The discrepancies become even larger in the presence 
of a constant magnetic field which magnifies the impor- 
tance of correlations. However, at high magnetic field 
QMC calculations are imprecise due to the fixed phase 
approximation and the exact diagonalization techniques 
require a number of configurations which increases expo- 
nentially with the number of particles, making the cal- 
culation unpractical for more than ten particles. An al- 
ternative method that has been shown to be quite accu- 
rate for cases in which HF yields broken symmetry so- 
lutions is that of symmetry restoration by projection^ 
although its complexity also increases very rapidly with 
the number of particles. The ground state properties 
of many electrons quantum dots can also be calculated 
with local-spin-density approximation (LSDA) and cur- 
rent density functional methods These approaches 
are based on energy functionals obtained through fits 
and interpolations of QMC results for the spin polar- 



ized and unpolarized two-dimensional electron gases that 
do not include a magnetic field and, moreover, provide 
a single particle spectrum whose physical interpretation 
is not clear. Nevertheless, LSDA has given good re- 
sults for quantum dots with low electronic densities, hav- 
ing Wigner-Seitz radius r s 1. These are the typi- 
cal densities of quantum dots made by etching 3 or by 
lithography 1 usually over a GaAs substrate. Another 
class of quantum dots are the self-assembled onesjSL 2 ^ 
with high electronic densities (r s < 0.5) and great prac- 
tical interest due to technological applications such as, 
e.g., in semiconductor lasers^ They are usually formed 
over InAs substrates, have small diameters and contain 
strongly confined electronic states. For systems so dis- 
tant from the homogenous conditions one should not ex- 
pect the LSDA to give accurate predictions and, indeed, 
most calculations on InAs dots have been made using 
non-DFT-based schemesi 30 i 31 i 32 

In this paper we employ the Brueckner-Hartree-Fock 
(BHF) metho d 33 ! 34 to describe electrons in quantum dots 
because it yields good ground state energies at a compu- 
tational cost that grows in a relatively modest way with 
the number of particles, as compared to CI and QMC (in 
its most accurate versions such as diffusion QMC). This 
method has been extensively applied in the past to study 
nuclei 3 ^ and metallic clusters^ Also, a related scheme 
based on Bethe-Goldstone equations has been used to 
study small excitonic complexes in quantum dots^i The 
underlying assumption of BHF theory is a description in 
terms of independent pairs of particles. For each pair, 
interaction affects the way in which the two particles 
scatter, while the rest of particles are considered non- 
interacting and their influence is only through the Pauli 
exclusion principle. It is well known that the exclusion 
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principle imposed by the Fermi sea induces a modifica- 
tion of the pair wave function at short interparticlc dis- 
tances, inducing the so-called "wound" in the wave func- 
tion. Therefore, BHF theory describes short range sta- 
tistical correlations and will miss, by construction, long 
range correlations associated with collective motions i^iM 
Here we compare BHF energies for electrons in quan- 
tum dots with the available results from methods at- 
tempting a direct solution of the many body Schrodinger 
equation like CI diagonalization and QMC, as well as 
with approximate models such as LSDA and HF. We find 
that BHF is always more accurate than HF. It recovers 
the exact result in the strong-confinement limit and at 
high magnetic fields, when the maximum-density-droplct 
(MDD) is formed, it yields more accurate energies than 
the LSDA ones. Furthermore, in the case of electronic 
densities typical of self- assembled quantum dots, BHF 
also gives much better results than LSDA even in the 
absence of magnetic fields. 



II. THE BHF SCHEME 

This section recalls some essential points of the BHF 
theory deriving several relevant equations and quoting 
others without proof. The restriction to circular sym- 
metry solutions and its practical implications are also 
discussed. 



The residual interaction V res fulfills 

H = HhF + V res , 
N N 

Vres = ^ Vi i ~y^, U ( i ) ~ C ■ 

i<j i 



(4) 
(5) 



HF theory yields the following general matrix elements 



{HF\V res \HF) = 
(HF\V r 
(HF\Vres\mni 



„— i 



'res\mi L ) = 



r 1 } * o , 



(6) 
(7) 
(8) 



where we have used the standard notation for particle- 
hole (ph) excitations, i.e., indexes i, j (m, n) refer to or- 
bitals below (above) the Fermi energy and \rni~ 1 } is the 
Slater determinant obtained promoting one electron from 
orbital i to orbital m in \HF). Note that only 2p-2h ex- 
citations yield non vanishing transition matrix elements 
since the two-body nature of V res ensures that matrix 
elements between determinants differing in more than 
two orbitals will again vanish. Another immediate con- 
sequence from Eqs. ©-© is that E HF = (HF\H\HF). 
We can write 



\*) = \HF) + J2*n\$n) ■ 
From Eqs. 101,0),© one easily finds that 



(9) 



A. Theory 

Let us consider a system of N electrons with Hamilto- 
nian 



jV N 
i i<j 



(1) 



where ho = U + v e xt{fi), ti is the single particle kinetic 
energy which can include the effect of the magnetic field, 
v ex t the confining potential and = e /e|rj — Tj\, with e 
the dielectric constant. The ground state |<J/) is solution 
of the many-body Schroedinger equation 



(H-E)\i£) = 



(2) 



where E is the ground state energy. Let us now con- 
sider an independent particle model (the HF model in 
the following) where the eigenstates are Slater de- 
terminants solutions of the equation 



(H HF -W n )\$ n ) = 0, 



(3) 



where Hhf = C + J^f ' {ha{i) + Ui), with Ui the HF poten- 
tial and C a constant. The HF ground state determinant 
and energy are given by |$o) = \HF) and Wo = Ehf, 
respectively. 



{Hhf — E) |HF) + ^ a „|$„) +V res |*) = 0. (10) 

V n^O J 

Multiplying Eq. JTJJ by (HF\ on the left one gets 

E = E HF + {HF\V res \$) . (11) 
If multiplying by ($ n | one finds a n — ^gj^™^ and 



E-Wn 



hence the following implicit equation is obtained 



\V) = \HF) + J2 



E-W n 



(12) 



This equation can be solved by iteration taking as start- 
ing energy the HF one: 



^) = \HF) + Y,-f 



QJVrrjHF) 



yielding for the energy 



Ehf - W n 



E — E HF + 22 



\(<S> n \Vres\HF)\ 2 

E HF - W n 



(13) 



(14) 



At the first order in V res this equation gives a result for 
the energy which coincides with the one of first order 
perturbation theory, summing all the orders we get the 
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correlation energy in the ladder approximation. This is 
more clear defining the G-matrix by the relation 



G\HF) = V res \V) 



(15) 



We then get the Bethe-Goldstone implicit equation for 
G: 



G = V„ 



^vJ^Ag , 



n#0 



E-W n 



(16) 



and from Eq. (|llfl 

E = E HF + (HF\G\HF) 
= E HF 



£ (HF\V na \$ n )($ n \G\HF) ( L _ ( 



E- W n 



Only 2p-2h determinants yield a non- vanishing contribu- 
tion to the sum of Eq. I|17fl . which can thus be reduced to 
a sum of two-body matrix elements. Assuming E = E^p 
on the right-hand-side, as in the ladder approximation, 
one has 



(ij\v\mn) 



e = e hf+ i^ _ 



tjmn 



((mn\g\ij) - (mn\g\ji)) , (18) 

where the e Q are the HF single particle energies and we 
have associated the G matrix with an effective two-body 
interaction g. 

In order to have a practical computational scheme it 
remains now to specify the two-body matrix elements 
of g in Eq. (JTSJ. This is accomplished within the BHF 
independent-pair model^ where the off-diagonal matrix 
elements arc found from 



[mn\g\i]) 



mn\v\ij) 

mn\v\pq)(pq\g\ij) 



E 



(19) 



The ground state energy of Eq. 1)18(1 with the matrix el- 
ements obtained from Eq. 1|19[) is the BHF energy which 
sums all the ladder diagrams corresponding to the iter- 
ated solutions of Eq. (|14fl . 



B. Circular symmetry restriction 

In this work we restrict to circular symmetry cases, 
where the HF orbitals can be factorized as 



(roji) = R nimi (r)—=x Pi (^) , 



2vr 



(20) 



where rii — 0,1,..., m, = 0, ±1, and \ii = ±1/2 
are the principal, angular momentum, and spin quantum 
numbers, respectively. In this situation the angular and 
spin parts of the two-body matrix elements yield selection 



rules on the corresponding quantum numbers and the 
matrix elements reduce to 



(ab\v\cd) = S ma+mb . mc+md 5^ c <J MI)W 



I r (R 



Rn 



3 5 Rn 



• ^n d m d ) i (21) 



where I r is a radial integral that we compute numerically. 
Note that through Eq. (|19|) the same angular momentum 
selection rules apply to (ab\g\cd) and that both matrix 
elements are real. 

The two-body matrix elements of g required for the 
evaluation of the total energy, Eq. (|18fl . are found by 
solving Eq. i|19|) as a linear system for the unknowns 
(mn\g\ij). For each pair ij we have an independent lin- 
ear system and the above mentioned selection rules are 
very important since they allow a big reduction in the 
number of effectively coupled equations. Since the space 
of particle states must be truncated, the convergence of 
the calculation with the number of empty HF states has 
to be controlled. Another check of the numerical accu- 
racy must be done regarding the number of radial points 
used in the evaluation of the integrals of Eq. I|21|l . 



III. RESULTS 

Table I compares the energies of B = ground states 
of iV-electron dots in BHF with the results of HF, 
LSDA (using the Tanatar-Ceperley parametrization for 
the correlation energy^) and diffusion QMCi^S The ex- 
ternal confinement is taken of parabolic type v ext (r) — 
miliar 2 /2, with m the electron effective mass. We refer all 
energies to the confinement energy Tiuiq and characterize 
the interaction strength by the repulsion-to-confinement 
ratio R, defined as 



R : 



V(e*o) 
hajQ 



(22) 



with £q indicating the oscillator length (Tiuo = h 2 /tti£q). 
The results in Table I correspond to R = 1.89. Tak- 
ing, for instance, typical GaAs values e—12A and m* = 
0.067m e the chosen R value would correspond to a con- 
finement energy of hu>o = 3.32 meV, which reproduces 
the experimental value of Ref. |3j. The number of elec- 
trons is varied from TV = 2 to 13. 

At B = the BHF energies obviously improve the HF 
ones, although they are still appreciably higher than the 
QMC and LSDA values. This is due to the fact that in 
BHF long-range collective correlations are missed. The 
importance of short range correlations is expected to in- 
crease as the system is more tightly confined, for a fixed 
strength of the Coulomb repulsion and, thus, a better 
performance of BHF is expected when increasing the con- 
finement strength. Indeed, this is shown to be the case 
in Fig. 1, where for N — 2 and 6 at B = we have 
varied the value of R. The lower panels of this Figure 
display the correlation energies, defined as usual by sub- 
tracting the HF value Ehf from the total energy E, i.e., 
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Ecorr = E — Ehf- Note that although the correlation 
energy is globally reduced when the ratio decreases BHF 
accounts for a larger part of it. Actually, for TV = 2 BHF 
accounts for 71% of the full (CI) correlation energy when 
R = 1 and 87% when R = 0.5. The corresponding figures 
for TV = 6 are 64% (R = 1) and 75% (R = 0.5). These 
numerical results are thus showing that in the limit of 
small R BHF converges to the exact correlation energy. 
It can be also seen from Fig. Q that at a given value 
of R, BHF correlations for the N = 2 dot are a some- 
what larger piece of the total correlation energy than for 
N = 6. We attribute this difference to a more impor- 
tant role of the collective effects leading to long-range 
correlations not included in BHF for the 6-electron dot. 

We focus next on the influence of a magnetic field and 
how the BHF energies are affected by it. In the sym- 
metric gauge, a magnetic field in the z direction (per- 
pendicular to the dot plane) induces a modification of 
the effective confinement from ujq to fl = ycJo + w%/4, 
with uj c the cyclotron frequency. For increasing magnetic 
fields one then expects short range correlations to be en- 
hanced due to the stronger confinement and, therefore, 
an improved performance of the BHF method. Figure El 
displays the evolution of the total energy with the ratio 
lu c /ujq for a fixed value of R = 1.5 and for a dot with 
2 electrons. Note that the correlation energy is about 
an order of magnitude higher for the singlet than for 
the triplet^ As expected, for increasing values of lo c /ujq 
BHF is accounting for a higher part of the correlation en- 
ergy, although the increase is rather moderate. For the 
singlet state BHF correlations range from 73% to 81% of 
the total correlation energy when uj c goes from to 5loq. 
The evolution is even flatter for the triplet, where BHF 
correlations remain at ss 75% for all values of oj c . Note 
also that due to the sizeable energy correction for the 
singlet, the singlet-triplet transition point is remarkably 
improved in BHF with respect to HF. 

Figure |21 shows the evolution with magnetic field of 
the results for a 6-electron dot with a fixed interaction- 
to-confinement ratio of R = 1.89. As compared to the 
N = 2 case, this dot shows a much richer phase dia- 
gram, with large variations in ground state angular mo- 
mentum and spin when increasing the magnetic field. 
Most remarkable is the comparison of BHF and LSDA 
energies: while LSDA is clearly superior to BHF at low 
fields the situation is reversed when entering the fully 
polarized phase corresponding to the maximum-density- 
droplet (MDD). In the MDD region the LSDA energy is 
actually higher than the HF one and only BHF is able to 
provide an energy correction in the right direction with 
respect to HF; total BHF energies being approximately 
halfway of CI and HF. 

All the BHF results shown above have been obtained 
using a large enough space of empty HF states, always 
checking that for the given accuracy convergence in Eq. 
(|18|l has been achieved. In practice, we include the low- 
est N p particle states and repeat the calculation increas- 
ing this number. Figure 0] shows the evolution of the 



BHF energy with N p , on scale proportional to l/N p , for 
three selected cases: the two upper panels correspond to 
six and nine electrons with the MDD configuration in 
strong magnetic field and a moderate confinement, while 
the lower panel shows the case of four electrons in very 
strong confinement and zero field. The results have been 
fitted with a polynomial including powers up to 
In the chosen scale, the N p — > oo limit is given by the 
intersection of the polynomial fit with the left vertical 
axis. It can be seen from Fig. 0] that the convergence 
is somewhat faster for the moderate confinement cases. 
However, in both examples the evolution with space di- 
mension is quite smooth, indicating that the correlation 
energy builds up by gathering contributions from many 
states. The N = 9 results of Fig. 01 provide additional 
support to our preceding conclusion from Fig.[3]that BHF 
performs much better than LSDA in the MDD region. 



As already mentioned, for strong confinements the 
BHF energies are close to the exact values. Indeed, the 
N — 4 results of Fig.0]are very illustrative in this respect 
since the extrapolated BHF energy essentially coincides 
with the CI result. This strongly confined system mimics 
a self-assembled InAs dot with Tiojq = 50 meV. Follow- 
ing Ref. |32Jwe take for this material to* = 0.024m. e and 
e = 15.15 giving, for N — 4, a small Wigner-Seitz radius 
of r s ~ 0.12. 



To emphasize the possibility of calculating the ener- 
gies of larger systems in BHF theory we end this sec- 
tion by showing in Fig. the results for an N = 20 dot. 
For this number of electrons exact methods like CI or 
QMC become extremely demanding and we have not at- 
tempted to compare with them. The evolution of the 
BHF energy with the number of particle states resem- 
bles that of Fig. 0] although, as one could expect, larger 
values of N p need to be considered for a similar degree 
of convergence. The proximity of the extrapolated-BHF 
and LSDA energies for N = 20 is a bit surprising since, 
as shown in Fig. [3J for 6 electrons in the same confine- 
ment the difference is larger. BHF correlations are thus a 
more important contribution for N = 20 than for N = 6. 
This can be understood as a different degree of magic- 
ity for these two dots. Indeed, a highly magic system is 
characterized by a distribution of single-particle orbitals 
whose energies group in bunches corresponding to quasi- 
degenerate shells, with large energy gaps between the 
shells. One expects a quenching of independent-pair mo- 
tions in a highly magic system, with respect to collective 
motions, and, therefore, a relatively worse performance 
of the BHF theory for them. Note also that since the 20- 
electron dot has a higher density than the 6-electron dot, 
for the same confinement, a better performance of BHF 
theory in the former agrees with our preceding results 
regarding the high density limit. 
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IV. CONCLUSIONS 

We have explored the prediction of BHF theory for 
the correlation energy in two-dimensional parabolic dots. 
Rather than attempting systematic calculations we in- 
tend to point out characteristic trends when applying a 
well established many body approach such as the BHF 
theory to parabolic quantum dots. By comparing with 
exact calculations we have quantitatively discussed the 
relevance of BHF correlations as a function of the con- 
finement potential and applied magnetic field for several 
quantum dots. BHF theory converges to the exact corre- 
lation energy in the limit of strong confinement potential 
(small R parameter, high densities). Also relevant is the 
limit of strong magnetic fields where the MDD is formed. 
In these two regions the standard LSDA is shown to be 
grossly inadequate while BHF stands as a competitive 
method that sizeably improves on HF theory. In general, 
the BHF correlation energy gathers contributions from 



many empty HF states, as evidenced by the smooth con- 
vergence with the space dimension. 

Possible extensions of the calculations presented here 
can consider i) relaxing the circular symmetry constraint, 
and ii) including self-consistency in the BHF single par- 
ticle orbitals. We shall address the latter by finding the 
improved mean-field proposed by Bethe, Brandow, and 
Petschek22 that actually implies the solution of a double 
selfconsistency problem, on orbitals and effective inter- 
action. Work along these lines is now in progress. 



Acknowledgments 

This work was supported by projects PRIB-2004-9765 
and FIS2005-02796 (LIS.). A. E. acknowledges support 
from the Generalitat de Catalunya. E.L. has been sup- 
ported by DGU (Spain) grant SAB2004-0091. 



* Permanent address: Dipartimento di Fisica, Universita di 
Trento, and INFN sezione di Trento, 1-38050 Povo, Italy 

1 U. Meirav, M. A. Kastner, and S. J. Wind, Phys. Rev. 
Lett. 65, 771 (1990). 

2 R. C. Ashoori, H. L. Stormer, J. S. Weiner, L. N. Pfeiffer, 
S. J. Pearton, K. W. Baldwin, and K. W. West, Phys. Rev. 
Lett. 68, 3088 (1992). 

3 S. Tarucha, D. G. Austing, T. Honda, R.J. van der Hage, 
and L.P. Kouwenhoven, Phys. Rev. Lett. 77, 3613 (1996); 
Jpn. J. Appl. Phys. 36, 3917 (1997); S. Sasaki, D. G. Aust- 
ing, and S. Tarucha, Physica B 256, 157 (1998). 

4 D.G. Austing et al, Phys. Rev. B 60, 11514 (1999). 

5 A. Kumar, S. E. Laux and F. Stern, Phys. Rev. B 42, 5166 
(1990). 

6 M. Fujito, A. Natori, and H. Yasunaga, Phys. Rev. B 53, 
9952 (1996). 

7 H. M. Muller and S. Koonin, Phys. Rev. B 54, 14532 
(1996). 

8 C. Yannouleas and U. Landman, Phys. Rev. Lett. 82, 5325 
(1999). 

9 B. Reusch and H. Grabert, Phys. Rev. B 68, 045309 
(2003). 

10 LI. Serra, R. G. Nazmitdinov, and A. Puente, Phys. Rev. 
B 68, 035341 (2003); A. Puente, LI. Serra, R. G. Nazmit- 
dinov, Phys. Rev. B 69, 125315 (2004). 

11 M. Eto, Jpn. J. Appl. Phys. 36, 3924 (1997). 

12 P.A. Maksym and T. Chakraborty, Phys. Rev. Lett. 65, 
108 (1990); 

13 D. Pfannkuche, V. Gudmundsson, and P. A. Maksym, 
Phys. Rev. B 47, 2244 (1993). 

14 P. Hawrylak and D. Pfannkuche, Phys. Rev. Lett. 70, 485 
(1993). 

15 J. J. Palacios, L. Moreno, G. Chiappe, E. Louis, and C. 
Tejedor, Phys. Rev. B 50, 5760 (1994). 

16 T. Ezaki, N. Mori, and C. Hamaguchi, Phys. Rev. B 56, 
6428 (1997). 

17 A. Emperador, E. Lipparini and F. Pederiva, Phys. Rev. 
B 72, 033306 (2005). 

18 A. Harju, V.A. Sverdlov, R.M. Nieminen, and V. Halonen, 



Phys. Rev. B 59, 5622 (1999). 

19 J. Shumway, L. R. C. Fonseca, J. P. Leburton, R. M. Mar- 
tin, and D. M. Ceperley, preprint cond-mat/9912166 

20 F. Bolton, Phys. Rev. B 54, 4780 (1996). 

21 C. H. Mak, R. Egger, and H. Weber-Gottschick, Phys. Rev. 
Lett. 81, 4533 (1998); R. Egger, W. Husler, C. H. Mak, 
and H. Grabert, Phys. Rev. Lett. 82, 3320 (1999). 

22 F. Pederiva, C.J. Umrigar and E. Lipparini, Phys. Rev. B 
62, 8120 2000. 

23 C. Yannouleas and U. Landman, J. Phys.: Condens. Mat- 
ter 14, L591 (2002); Phys. Rev. B 69, 113306 (2004). 

24 M. Koskinen, M. Manninen, and S.M. Reimann, Phys. 
Rev. Lett. 79, 1389 (1997). 

25 K. Hirose and N. S. Wingreen, Phys. Rev. B 59, 4604 
(1999). 

26 M. Ferconi and G. Vignale, Phys. Rev. B 50, 14722 (1994); 
E. Lipparini, N. Barberan, M. Barranco, M. Pi, and LI. 
Serra, Phys. Rev. B 56, 12375 (1997); M. Pi, M. Barranco, 
A. Emperador, E. Lipparini, and LI. Serra, Phys. Rev. B 
57, 14783 (1998); O. Steffens, U. Rossler, and M. Suhrke, 
Europhys. Lett. 42, 529 (1998); 44, 222 (1998). 

27 B. T. Miller, W. Hansen, S. Manus, R. J. Luyken, A. Lorke, 
J. P. Kotthaus, S. Huant, G. Medeiros-Ribeiro, and P. M. 
Petroff, Phys. Rev. B 56, 6764 (1997). 

28 R. J. Warburton, B. T. Miller, C. S. Diirr, C. Bodefeld, K. 
Karrai, J. P. Kotthaus, G. Medeiros-Ribeiro, P. M. Petroff, 
and S. Huant, Phys. Rev. B 58, 16221 (1998). 

29 D. Bimberg, M. Grundmann, N. N. Ledentsov, Quan- 
tum dot hetero structures (John Wiley and Sons, England, 
1999), Chap. 8. 

30 B. Szafran, J. Adamowski, and S. Bednarek, Phys. Rev. B 
61, 1971 (2000). 

31 S.-J. Cheng, W. Sheng, and P. Hawrylak, Phys. Rev. B 
68, 235330 (2003). 

32 T. Brocke, M. T. Bootsmann, M. Tews, B. Wunsch, D. 
Pfannkuche, Ch. Heyn, W. Hansen, D. Heitmann, and C. 
Schiiller, Phys. Rev. Lett. 91, 257401 (2003). 

33 M. A. Preston, R. K. Badhuri, Structure of the nucleus 
(Addison- Wesley, USA, 1982), Chap. 8. 



6 



34 A. L. Fetter and J. D. Walecka, Quantum theory of many- 
particle systems (McGraw-Hill, USA, 1971). 

35 K. A. Brueckner, J. L. Gammel, and H. Weitzner, Phys. 
Rev. 110, 431 (1958); J. W. Negele, Phys. Rev. C 1, 1260 
(1970). 

36 E. Lipparini, LI. Serra and K. Takayanagi, Phys. Rev. B 
49, 16733 (1994). 

37 R. Perez and A. Gonzalez, Physica E 8, 91 (2000). 

38 B. Tanatar and D. M. Ceperley, Phys. Rev. B 39, 5005 
(1989). 

39 H. A. Bethe, B. H. Brandow, and A. G. Petschek, Phys. 
Rev. 129, 225 (1963). 



Tables 




TABLE I: Ground state energies for the dots with 2 < N < 
13 computed by HF, BHF, LSDA, QMC and CI methods. 
The energies E in units of the confinement energy Tiluq axe 
tabulated as E' = E/hcoo- A fixed value R = 1.89 of the 
interaction-to-confinement ratio has been used. 



N 




^BHF 


-E"LSDA 


^QMC 


E'ci 


2 


4.078 


3.832 


3.739 


3.650 


3.646 


3 


8.589 


8.289 


8.082 


7.979 


7.957 


4 


13.94 
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Figures 



FIG. 1: (color online) Upper panels: total energies in dif- 
ferent methods for 2 and 6-electron dots as a function of the 
interaction-to-confinement ratio (see text). Lower: Correla- 
tion energies within each model for the same two dots. 
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FIG. 2: (color online) Left (right) panels show total (corre- 
lation) energies for the models and states indicated by the 
corresponding labels. The results correspond to the N = 2 
dot in a magnetic field, shown as a function of the cyclotron 
frequency (in units of u>o). The interaction-to-confinement 
ratio R (see text) is chosen as R — 1.5. 
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FIG. 3: (color online) Evolution of the total energy of a 6- 
electron dot with the magnetic field. The different phases are 
indicated by the angular momentum labels (L Z ,S Z ). A fixed 
value of the parameter R = 1.89 has been used (see text). 
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FIG. 4: (color online) Evolution of the total BHF energy with 
the number N p of empty HF states included in the solution of 
Eq. 1191 . Each panel shows the results for a different quantum 
dot. The solid line is a cubic fit, in powers of 1/N P , allowing 
extrapolation to the N p — > oo limit. The LSDA, HF and CI 
energies for each case are also shown for comparison. 
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FIG. 5: (color online) Same as Fig. 0] for a dot with 20 elec- 
trons and the additional parameters given in the inset. 
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